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Abstract: We propose two experimental schemes that can produce 
an arbitrary photon-number entangled state (PNES) in a finite dimen- 
sion. This class of entangled states naturally includes non-Gaussian 
continuous-variable (CV) states that may provide some practical advantages 
over the Gaussian counterparts (two-mode squeezed states). We partic- 
ularly compare the entanglement characteristics of the Gaussian and the 
non-Gaussian states in view of the degree of entanglement and the Einstein- 
Podolsky-Rosen correlation, and further discuss their applications to the C V 
teleportation and the nonlocality test. The experimental imperfection due to 
the on-off photodetectors with nonideal efficiency is also considered in our 
analysis to show the feasibility of our schemes within existing technologies. 
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1. Introduction 

Ever since Einstein-Podolsky-Rosen (EPR)'s argument against quantum mechanics was put 
forward [1 1, quantum entanglement has been a topic of great interest from a fundamental point 
of view. It can play a crucial role in manifesting striking differences between quantum and clas- 
sical (e.g. local hidden-variable |2}-|4]) descriptions of nature. Furthermore, it has also drawn 
much attention of practical interest because quantum correlations can be employed to carry out 
information tasks to the extent far beyond their classical counterparts, e.g. quantum comput- 
ing (5) and teleportation (6). In a bipartite setting, the primitive entangled states for discrete 
variables are the so-called Bell states, the maximally entangled states of two qubits, e.g. singlet 



state. In the regime of continuous variables (CVs), the Bell state can be realized in the form of 
two-mode squeezed state (TMSS), which becomes maximally entangled in the limit of infinite 
squeezing. The TMSS has been mostly the target entangled resource to produce for various 
quantum information tasks Q like C V quantum teleportation jH ■ 

The TMSS belongs to the class of Gaussian states, which has been extensively studied both 
theoretically and experimentally for CV quantum informatics J9). On the other hand, a great 
deal of attention has also been directed to the non-Gaussian regime (e.g. state engineering ifTUI 
and characterization ifTTlfTSl '). as the non-Gaussian entangled states can provide some prac- 
tical merits lfT3Tl2Tl and even become an essential ingredient 12214251 for a number of quan- 
tum tasks. Furthermore, when the quantum information processing is performed under realistic 
conditions, the quantum correlation is inevitably degraded and it thus becomes an important 
question whether Gaussian or non-Gaussian entangled states can be more robust against deco- 
herence 1261271 . It was recently demonstrated that there exists a broad parameter space in which 
non-Gaussian entanglement can survive longer than Gaussian entanglement under noisy envi- 
ronments 1 28 30] or quantum-limited amplifier Q"T1 . With all these considered, it seems very 
desirable to have an experimental toolbox to generate a broad class of non-Gaussian entangled 
states in a controllable way. 

In this paper we consider a class of CV entangled states in the photon-number entangled form 
Ln=o Cn \n) a \n)b, where \n) denotes a Fock state basis. One particular example is the TMSS with 
the coefficients C„ = A"(l — A 2 ) 1 / 2 (A: squeezing parameter, N — > °°), which is the only Gaus- 
sian state among the photon-number entangled states. Another example is the pair-coherent 
state given by C„ ~ C"/ n ' E2 33] which can be useful for a number of applications includ- 
ing quantum teleportation 041 . quantum metrology 031 . and a Bell test Q~6l . In fact, a broad 
class of photon-number entangled states has been so far considered for the nonlocality test us- 
ing homodyne detections 0~7] - l39l . Here we propose two experimental schemes to generate a 
finite-dimensional PNES with arbitrary coefficients, Y^=oC n \n) a \n) y (X^ = o|C ; ,| 2 = 1), where 
the coefficients C„ can be controlled with beam splitting and squeezing parameters. Both of 
our proposed schemes make use of coherent superposition operations in single-photon interfer- 
ometic settings that erase the which-path information on the realized photonic operations. The 
first scheme employs the second-order superposition operation tad 1 ' + ra T a, which has been 
recently proposed and experimentally implemented in the context of proving bosonic commu- 
tation relation [a, a'] = 1 B014431 . together with two-mode squeezing operation. We note that 
the coherent operation taa' + ra^a was also discussed in the context of noiseless quantum am- 
plifier 1441 . On the other hand, the second scheme employs a sequence of nonlocal first-order 
coherent superpositions ta + rb\ Its single-mode version ta + ra' was recently proposed for a 
quantum state engineering 11451 and also shown to be useful to enhance two-mode entanglement 
properties 112011461 . 

We also address the usefulness of the finite-dimensional PNES for CV quantum teleportation 
J8) and nonlocality test 114711481 compared with the two-mode squeezed state. Furthermore, 
it was very recently shown that the photon-number entangled states in finite dimension, e.g. 
Co|0) ( ,|0)/, + Ci I can survive longer under noisy environments than the TMSS with the 

same degree of entanglement or energy lf28]|29l . Therefore, our proposed schemes can be a 
useful tool not only for CV quantum applications but also for fundamental tests of quantum 
physics. 

This paper is organized as follows. In Sec. 2, we first compare the entanglement properties of 
finite-dimensional PNES (non-Gaussian) and a two-mode squeezed state (Gaussian) in view of 
the degree of entanglement and the EPR correlation. In Sec. 3, we further investigate the use- 
fulness of the PNES for CV quantum teleportation and nonlocality test. Then, we propose two 
experimental schemes to generate an arbitrary PNES, Y^=oCn\ n )u\ n )b, m Sec. 4. We illustrate 




the feasibility of our schemes in Sec. 5 by investigating the generation of a PNES up to two- 
photon correlation, i.e. Cb|0) a |0)ft + Ci|l) a |l)j,+C2|2) a |2)j, considering realistic experimental 
conditions. In Sec. 6, our results are summarized. 

2. Entanglement and EPR correlation 

First, we briefly compare the TMSS and the PNES in terms of entanglement properties in or- 
der to identify the practical relevance of the PNES for CV quantum informatics. For a pure 
two-mode state | 1 P)ab, the degree of entanglement can be quantified by the von Neumann 
entropy E(pa) = — Tr^ [p^ log 2 Pa] for the reduced density operator p A = TibW^ab^Iab]- 
For the class of photon-number entangled states Ym=oCn\n)a\n)b> the von Neumann entropy 
becomes maximal when all the coefficients C„ are identical. Thus, for the case of TMSS 
with C„ = A"(l — A 2 ) 1 / 2 , the state can have an infinite degree of entanglement with infinite 
squeezing, i.e. A = 1, which is practically impossible to achieve. On the other hand, the finite- 
dimensional PNES can match or even surpass a finitely-squeezed TMSS as an entangled re- 
source. In Fig. 1(a), we plot the degree of entanglement for the TMSS (blue solid) as a function 
of the squeezing parameter s = tanfr'A. This is compared with the maximal possible en- 
tanglement for the PNES with equal coefficients {C\ = C2 = ■ ■ ■ = Cn) of dimensions N = 1 
(red dotted), 2 (red dashed), and 10 (red dot-dashed). The degrees of entanglement for the 
PNESs are given by 1, 1.585 and 3.459, respectively. To achieve such degrees of entangle- 
ment, the squeezing of the TMSS should be s = 0.5185 (4.506 dB), s = 0.7335 (6.374 dB) 
and s = 1 .391 (12.09 dB), respectively. In the pulsed-regime generation of squeezed states, the 
level of squeezing currently available from an optical parametric amplifier is s = 0.403 (3.5 dB) 
J49]|50l so that the PNES with N — I can already surpass the entanglement of the TMSS. 

We also look into another entanglement property, the EPR correlation, which is the total 
variance of a pair of EPR-like operators, EPR = A 2 (xa — xb) + A 2 (pA + Pb)- Here Xj = -j- (dj + 

at) and pj — -^(fl/ — flj) (j = A,B) are the quadrature amplitudes of the field that can be 
measured in homodyne detection. The value of EPR below 2 represents the quantum correlation 
between the quadrature amplitudes of two modes. In Fig. 1(b), the EPR correlations of the 
PNESs for the dimensions N = 1, 2, and 10 are 1.172, 0.8315, and 0.2516, respectively. The 
corresponding levels of squeezing for the TMSS are given by s = 0.2674 (2.324 dB), s = 
0.4388 (3.813 dB) and s = 1.037 (9.008 dB). Thus the PNES with N >2 can surpass the 
currently available TMSS (s = 0.403) in view of the EPR correlation. 



3. Applications: CV teleportation and nonlocality test 

In this section, we further investigate the usefulness of a finite-dimensional PNES particularly 
for continuous variable (CV) teleportation and nonlocality test. For this purpose, we evaluate 
the quality of CV teleportation by the average fidelity between an unknown input state and the 
teleported state [ 8] , and investigate the nonlocality test by Banaszek and Wddkiewicz based on 
the phase-space distribution functions 147 48l. 

(i) The teleportation fidelity in the Braunstein-Kimble (BK) scheme |8) can be evaluated in 
terms of the characteristic functions of an input state and its teleported state as 

F = - />AC out (A)C m (-A), (1) 
n J 

where C ou t(A) = Ci n (A)C£(A*, A) l53l . Here Ce{X*,X) is the characteristic function of a two- 
mode entangled state. We consider the finite-dimensional PNES L^=o \ n )a I") b f° r the dimen- 
sions N = 1,2 and 3 as an entangled resource. For instance, the characteristic function of the 
PNES for N = 2 is given by 

C £ (A 2 ,A 3 ) = e -(^l 2 +N 2 )/ 2 [|Cb| 2 + |d| 2 (l - |A 2 | 2 )(1 - |A 3 | 2 ) 
+ M(2 -4\X 2 \ 2 + |A 2 | 4 )(2 -4|A 3 | 2 + |A 3 | 4 ) 

„ 1*1* , r< r«f\ 1 i 0)^2 o *2o *2 i CqC 2 ,2i2 
+ Cq(_1 A 2 A3 + A2A3 i — A2 A3 -\ — A2A3 

+ i(CtC 2 A 2 *A3* + C 1 C 2 *A 2 A3)(|A 2 | 2 -2)(|A3| 2 -2)], (2) 

where |C | 2 + \Q | 2 + |C 2 | 2 = 1. For the case of teleporting an arbitrary coherent-state input, we 
find, by optimizing the fidelity (1) using Eq. (2), that the average fidelity can be achieved up to 
F = 0.7334 at the choice of Co « 0.765, C\ « 0.535 and C 2 « 0.359. 

In Fig. 2(a), we compare the teleportation fidelity achieved via the PNES and the fidelity via 
the TMSS. The optimal fidelity via each PNES at N = 1, 2 and 3 corresponds to the fidelity 
via the TMSS with the squeezing parameters s = 0.320 (2.776 dB), s = 0.506 (4.397 dB), and 
s = 0.638 (5.548 dB), respectively. Thus, the PNES at N = 2 can surpass the fidelity via the 
TMSS with the currently available squeezing in the pulsed regime, i.e., s w 0.403 (3.5 dB) l53l . 
As we show in Section 4, our proposed schemes do not require a high-level of squeezing to 
produce the optimal PNES (case of N = 2) for CV teleportation. 

(ii) We next consider the nonlocality test by Banaszek and Wodkiewicz that is addressed in 
phase space using the two-mode Wigner function [47 48 1. This Bell inequality is given by 

2 

Bw = ^-\W(a,P) + W(a,p') + W(a',P)-W(a',p')\<2, (3) 

where W(a, J3) is the two-mode Wigner function. We find that for the PNES, Y^,=(,C„\n) u \n) b at 
N = 2, the Bell inequality can be violated up to Z?bw = 2.32088 with the coefficients Co ~ 0.589, 
C\ w 0.700 and C 2 w 0.404. In Fig. 2(b), we see that this degree of nonlocality almost reaches 
the level for the TMSS with infinite squeezing l54l . Furthermore, we also see that the Bell 
violation by the PNES at N = 3 surpasses the value Z?bw of the TMSS in the entire region of 
squeezing. We can again achieve such degree of Bell violation using the weak squeezing regime 
in our schemes (Sec. 4). 

4. Experimental schemes 

We now propose two optical schemes to generate an arbitrary PNES, Y^!=o^"\ n )a\ n )b- O ne 
scheme is based on a second-order coherent superposition operation with two-mode squeezing 




Fig. 2. (a) Average fidelity in teleporting a coherent state and (b) Bell parameter Bbw 
as a function of the squeezing parameter s for the \TMSS) (blue solid) and the PNES 
Ln=o C "\ n )a\ n )b at N = 1 (red dotted), N = 2 (red dashed) and N = 3 (red dot-dashed). 
The coefficients of the PNESs are optimized for each N, 



operations, and the other on a sequence of first-order coherent superposition operations. 

(i) We first consider the operation tad 1 ' + rd r d acting on a single-mode a, which is the coher- 
ent superposition of two product operations — photon addition followed by subtraction (ctfft) 
and photon subtraction followed by addition (d^a). This coherent operation was experimen- 
tally implemented to prove the bosonic commutation relation [d,d'] = 1 [40,41]. While such 
a commutator is addressed as an equal superposition of the two product operations, dd r — d^d, 
we adopt an arbitrarily weighted superposition of the two operations, i.e., tad 1 ' + rcftd. In par- 
ticular, we show that the single-mode operation tad ' + rd^d together with two-mode squeezing 
operations S a b(%) = ex P(~ ^d 1 }? 1 + ^*bd) can constitute an essential building block to generate 
an arbitrary PNES. 

Suppose that a two-mode squeezing S a b(£)> the coherent operation faa^ + ra'd and the in- 
verse squeezing SL,(<jj) are sequentially applied to an input state. That is, we apply a sequence 
of operations defined by 

O" =^lb(^n)(t n dd 1 ' + r n d A a)S ab (E, n ) 

=A„ + (t„ + r„)(fl T acosh 2 s n +^Z?sinh 2 .s', 7 ) 
— (t„ + r„) coshs„ sinhs„ [exp(— iq> n )db + exp(/<p„)a^] , (4) 

where 

A„ = t„ cosh s„ + r n sinh" s„ , (5) 

with |f„| 2 + \r n \ 2 = 1. In the above Eq. (4), the identity Sl b (%)dSai,(£) = acoshi — 
^exp(;<p)sinh.y and §l b (£,)d^S a t,(£,) = a^coshs — /?exp(— z<p) sinhs are used, where £, = 
sexp(iq>) I531 . 

When a vacuum state |0) a |0)& is injected as an input, Oi yields a superposition of number 
states as Oi|0) a |0)f> = cosh 2 si[(fi +ntanh 2 si)|0)„|0)fo-exp(;'<Pi)(fi + n)tanhsi|l) a |l)/,]. In 
principle, a succession of 0„ applied on the vacuum states, n^=i O«|0) a |0)j, can yield an ar- 
bitrary superposition state Ym=o^"\ n )o\ n )b by properly choosing the parameters s n , r n , t n and 
(p„. For instance, the state 0\ |0) ( ,|0)/, ~ Cb|0)„|0)/, +Ci|l) a |l}i, can have a larger proportion 
°f |l)a|l)fc> i- e - | Co | < |Ci|, under the condition ritanhsj > t\. For comparison, if one instead 
applies the original quantum scissor scheme on the TMSS that projects an input onto the sub- 
space spanned by |0) and |1) [56], the output state becomes ~ |0)„|0)i + A|l) a |l)j,. That is, 




Fig. 3. (a) Experimental scheme to implement the operation S^ b (^)(taa [ + ra' a)S a t,(^) on 
an arbitrary state. BS1, BS2, and BS3 are beam splitters with transmissivities T\, Ti and t„, 
respectively. PDO, PD1 and PD2: photo detectors. The operation is successfully achieved 
under the detection of a single photon at only one of two detectors PD1 and PD2, with 
PDO clicked, (b) For a vacuum input state, the sequence of operations O n can yield a finite 
dimensional PNES, Y^ n= QC n \ri) a \n)i,. 



the vacuum state |0) a |0)f, is always more weighted than the single-photon state On 
the other hand, a generalized scissor scheme proposed for the noiseless quantum amplifier ll57l 
can be used to control each coefficient arbitrarily in the |0)-| 1) subspace, which will be briefly 
discussed in Sec. 5. 

The elementary operation 0„ can be experimentally realized as depicted in Fig. 3(a). An input 
state \y)ab is first injected into a nondegenerate parametric amplifier (NDPA) with coupling 
parameter <!;„, and then into the beam splitter BS1 (transmittance: T\ w 1) with the other input 
mode c in a vacuum. This can be described by 

&acSabgn)\V)ab\0)c « (1 - ^-ac f )S ab ^ n )\w)ab\0)c (6) 

Next the mode a is further injected into another NDPA with small coupling s <C 1 and the output 
is kept only under the condition of single-photon detection at PDO. That is, 

e(i\S ae (l - ^-ac^S ah ^ n )\w)ab\0) c \0) e « ~sa\l - ^K t )$*(6,)|Y'U|0)c, (7) 
which is then injected into the BS2 (T2 ~ 1), 

a* 

(s)B ad aHl - Tf-^)Sab{in)W)ab\0)cd 

« (s){\ - %<P)a\l - ^-atf)S ab £ n )\y) ab \0) cd , (8) 



where |0)„/ = |0) c |0)</. The next beam splitter BS3 making the transformations c — > t„c + r„d 
and d — > t„d — r n c gives 

\S\ W) ) = (s)[l -^fllf„^-r„c t )]a t [l--lfl(f, I c + + r„i + )]5 flfc (^)|^) afe |0)« ; . (9) 
i 2 

On detecting a single photon at PD1 (PD2) and no photon at PD2 (PD1), the state is projected to 

j. j_ /\ ft* ft* ft* ft* 

{tda* + rd" d)S a b{^n)W)ab^ where t ~ j^sf,, and r ~ ^ir,, (f ~ —-^sr n and r ~ j^^„)- Finally, 
the NDPA with the coupling parameter — % n yields |v) ou t ~ ^((^(faa 1 + ra t a)S„£ (<!;„) | v) a i,, 
with the identity S\ b {t, n ) = S ab (-£, n ). 

(ii) Second, we show that the sequence of two first-order coherent superposition operations, 
(t 2n a + r 2n P)(t 2 n _iZ? + r 2n _ia' ), can also yield an operation similar to 0„ in Eq. (4). A similar 
type of coherent operation was previously investigated in a form acting on a single-mode, td + 
ra\ which is the superposition of photon subtraction and addition l45l . Here we consider a 
nonlocal coherent superposition acting on two modes, td + rb' (tb + ra*). 

We define an operator 

On = (t2n« + ) (fen- 1 b + r 2n - lfl f ) 

=t2n- 1 hndb + r 2n - ir 2n fl t & + r 2n - 1 hnda* + t 2n -i nJ^b, ( 1 0) 

where |t 2n -i | 2 + |r 2n _i | 2 = 1 and |t 2n | 2 + |r 2n | 2 = 1. Given a vacuum state as an input, O l yields a 
superposition of number states as 0\ |0) a |0)& = ri (t 2 |0) a |0)fc + r 2 |l) a |l)j,). Furthermore, a suc- 
cession of O n , i.e., n^=i O n \0) a \0)b, can yield any desired superposition state Y*n={)C n \n) a \n) b 
by properly choosing the parameters r 2n _i, t 2n _i, r 2n and t 2n . Here the coefficients can be read- 
ily controlled only by the beam-splitter parameters as shown below. 

The operation O n can be implemented as depicted in Fig. 4. First, an arbitrary two-mode state 
| y/} a b is injected into an NDPA with small coupling si <C 1 and a BS1 with high transmissivity 
Ti«l, with mode a (b) into NDPA (BS1). The other input modes to the NDPAs and the BSs 
are all in vacuum states. Then, the BS3 (transmissivity: t 2n _i) yielding the transformations 
c T — > t 2n _ic^ +r 2n _iii^ and d' — s- t 2n _id' — r 2n _ic T gives the output 

R* 

[1 -^^(t2n-l^ t -r 2 n-lC + )][l-Sifl t (t2n-lC + +r2n-li + )]|r)«fc|0) £ - ( /- (H) 
M 

With the detection of single photon at PD1 (PD2) and no photon at PD2 (PD1), we see from 
Eq. (11) that the state is projected to \<&)ab = fe n -i^ + r ?n-i^ + )l W)ab> where t 2n _ 1 ~ — ^fen-i 

R* ' 

( Yj-r 2n _i) andr 2n _! ~ — sir 2n _i (— Sit 2n _i). Next, the output state is further injected into another 
NDPA with small coupling s 2 <C 1 and a BS2 with high transmissivity T 2 rj 1, with mode a (b) 
into BS2 (NDPA). Finally, a beam splitter BS4 (transmissivity: t 2n ) yielding the transformations 

e f -> t 2n e f +r 2n / t and / f -> t 2n / + -r 2n e f gives 

I^^EEll-^fl^+r.n/^lIl-S^fenf-r^)]^)^^),/. (12) 

Once again, with the detection of single photon at PD3 (PD4) and no photon at PD4 (PD3), 
we see from Eq. (12) that the state is projected to (t 2n a +r 2n £> T ) (t 2n _i^ + r 2n -i^ T )l V)ab, where 

1 R* R* / 

l 2n % r 2n (- xf l 2n) and r 2n s 2 t 2n (s 2 r 2n ). 

(Hi) In Sec. 3, we have seen that the optimal PNES T^ =0 C n \n) a \n) b with N — 2 for CV 
teleportation has the coefficients Co ~ 0.765, C\ m 0.535, and C 2 ~ 0.359. Under our first 
scheme, these coefficients can be obtained using the experimental parameters, e.g. s\ = ,s 2 = 0.1, 




Fig. 4. Experimental scheme to implement the operation (t 2n <3 + r^fo'' )(t2n-l^ + r 2n-l«' ) 
on an input state \ \ff) a b- BS1, BS2, BS3 and BS4 are beam splitters with transmissivities Tj, 
T 2 , i2 n -i, an£ i fem respectively. PD1, PD2, PD3 and PD4: photo detectors. The operation 
is successfully achieved under the detection of a single photon at only one of two detectors 
PD1 and PD2 and the detection of a single-photon at only one of two detectors PD3 and 
PD4. 

01 = 0, 02 = n, n « 0.4589, and r 2 « 0.9984, with t x = (1 - r\ ) 1 / 2 and t 2 = -(1 - r\) l l 2 . 
On the other hand, under the second scheme, the same coefficients can be obtained using the 
parameters r 2 « 0.3, r 3 sa 0.3863 and r 4 w 0.6193, with t 2 = -(1 - r^) 1 / 2 , t 3 = (1 -r 2 ,) 1 / 2 , and 
t4 = (1 -r|) 1/2 . Note that in our first scheme, we generate such a PNES using the NDPAs in 
the weak squeezing regime, s — 0.1 (0.869 dB). Furthermore, in the second scheme, we can 
obtain the same PNES only by adjusting the beam-splitter parameters, therefore, a high-level 
of squeezing is not necessary in our schemes. This is also true for the case of nonlocality 
test shown in Sec. 3. Under the first scheme, the optimized coefficients of the PNES for the 
nonlocality test are obtained using the parameters s\ = s 2 = 0.1, (pi = 0, </> 2 = r\ ~ 0.38 
and r 2 « 0.999, with t\ = (1 — r 2 ) 1 / 2 and f 2 = (1 — r 2 ) 1 / 2 . Under the second scheme, they 
are obtained using the parameters r 2 w 0.3, r 3 w 0.391 and r 4 rj 0.221, with t 2 = (1 — rj) 1 / 2 , 
t 3 = (l-r 2 ) 1 /2 i andt 4 = (l-r 2 ) 1 /2. 

5. Experimental feasibility 

In order to further illustrate the feasibility of our proposed schemes, we address realistic exper- 
imental conditions in producing a PNES up to one-photon correlation, Co |0)„ |0) i, + Ci |l) a |l)j, 
or two-photon correlation, Co|0) a |0)& + Ci|l)a|l)f> + C 2 |2) a |2}/,, as examples. In the two 
schemes of Figs. 3 and 4, we particularly consider each photodetector as an on-off detector 
that only distingushes two events, detection and non-detection, with efficiency 77 . The photode- 
tection can then be characterized by a positive operator-valued measure (POVM) [58 - 60 1 with 
two components Flo = — (n\ (no click) and fli =7 — fto (click). 

In the first scheme, an arbitrary input state goes through a sequence of operations — a two- 
mode squeezing, a second-order superposition tact' + rata heralded by nonideal on-off de- 
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Fig. 5. Fidelity between the ideal state Co|0} a |0}i, + Ci|l) a ]l)i> and the output state p ou t 
obtained by applying § ab (^)(tad' + ra' d)S a t,(t;) (blue circle) or (t2fl + r2^ ){hb + r l«' ) 
(red square), using on-off detectors with efficiency rj to the input state p m = |0) a |0)b as a 
function of | Cq | 2 for r\ = 0.66. Black triangle represents the output fidelity using the scissor 
scheme of 1571 . with the input two-mode squeezed state (s = 0.1) and the on-off detectors 
(7] = 0.66). 



tectors, and the inverse two-mode squeezing. This sequence yields an output state 

Pout— -7. T-77 7. r+T; (i-J) 

Twjngnfnf^pi,,^] 

where p in = |v)(vU ® \ty{0\cde and (7i = sl b B cd B ad SaeB ac S ab . We can evaluate the perfor- 
mance of our scheme by investigating the fidelity Fn = (V'a'IPou? I Wn) between the ideal target 
state, I y/f/) = £^ =0 C„|n) a \n) t, (N = 1,2) and the corresponding output state, pj^t (N — 1,2). In 
Fig. 5, we first show the case of PNES up to one-photon correlation, where the fidelity F\ (blue 
dot) is plotted as a function of |Co| 2 with the detector efficiency 77 = 0.66. With the parameters 
si =0.1 and 7^ = T 2 2 = 0.99 in Fig. 3, we see that a high fidelity above 0.996 is achieved in 
the whole range of |Co| 2 , with the detector efficiency 77 = 0.66 currently available i6TH63l . For 
comparison, we plot the fidelity of the output state using the generalied scissor scheme of ll57l 
with the input two-mode squeezed state (s = 0.1) and the on-off detectors (rj = 0.66). We see 
that our schemes yield a slightly better fidelity than the scissor scheme. 

In Fig. 6, we further investigate the fidelity F2 for the case of PNES up to two-photon cor- 
relation as a function of |Ci| 2 and | C2 1 2 - With the same parameters (tj = 0.66, s\ = 0.1 and 
r 2 = T 2 2 = 0.99) used in Fig. 5, the fidelity is achieved at least above 0.941 in the whole range 
of I Ci 1 2 and IC2I 2 . For both of the cases, the fidelity slightly increases with the vacuum-state 
probability |Co| 2 , as the weak coupling (si = 0.1) of the NDPA makes a low photon-number 
state better controlled. 

In the second scheme, two first-order superposition operations (t2« + 1"2^^) (ti & + ri ) her- 
alded by nonideal on-off detectors are sequentially applied to an arbitrary input state. This 
yields an output state 

n , Tr^ingnfngnf^p;^] (u) 




Fig. 6. Fidelity between the ideal state Co|0} a |0}fe + Ci|l) a |l}i>+C2|2) a |2)j and the output 
state Pout obtained by applying twice (a) S f ab (^ 2 )(t2aa' f + r 2 a' f a)S ab {£, 2 )sl b {^i){tida^ + 
ria f ci)S a t,{£,i) or (b) (t^ + ^S^^S + ^a^^a + ^B^X^S + r!^"), using on-off detectors 
with efficiency rj to the input state p; n = |0} a |0)j as a function of |Ci| 2 and \C 2 \ 2 for 
rj =0.66. 



where p- n = | v)(vU ® \0){0\cdef and U 2 = B ef B ae S hf B cc iB bd S ac . We investigate the fidelity 
f n = (VzvlPout I Vat) between the ideal target state, |^) = I^=oC„\n) a \n) b (N = 1,2), and the 
corresponding output state, p [^ (N = 1,2). In Fig. 5, we plot the fidelity Fj (red square) as 
a function of |Co| 2 with the detector efficiency rj = 0.66. With the parameters Si = S2 = 0.1 
and T 2 = l\ = 0.99 in Fig. 4, we find that a high fidelity above 0.993 is achieved in the whole 
range of |Co| 2 . In Fig. 6, we investigate the fidelity F 2 as a function of |Ci| 2 and |C 2 | 2 . With 
the same parameters (rj = 0.66, si = s 2 = 0.1 and T 2 = T 2 = 0.99) used in Fig. 5, the fidelity 
is achieved at least above 0.949 in the whole range of |Ci| 2 and |C 2 | 2 . Therefore, both of our 
schemes seem to make an output state at a very high fidelity even with nonideal on-off detectors 
used for heralding the conditional generation of the PNES. 

We have also calculated the success probability numerically for the output states under each 
scheme. For the state Co|0) fl |0)& +C\ | the first scheme, with the condition of s\ =0.1, 

s = Q.\ and T 2 = T 2 2 = 0.99, yields the success rate in the range of 2.4 x 10~ 6 (|C | 2 = 1 /2) 
to 10~ 4 , which increases with the coefficient |Cb|. On the other hand, the second scheme, with 
the condition of S] = s 2 = 0.1 and T 2 = T 2 , = 0.99, yields the success rate ~ 10 -4 . The success 
probability can of course be made larger by using higher-squeezing NDPAs in each scheme at 
the expense of output fidelity to some extent. 

Other than nonideal detector efficiency, dark counts might potentially degrade the output fi- 
delity. However, a recent experiment reported that a coincidence detection scheme recording 
only the synchronized events of laser pulse and a detector click in the pulsed regime can sig- 
nificantly eliminate dark count events l64l . Another experimental imperfection may also arise 
from the error in the transmissivity of beam splitter under our proposed schemes. In Fig. 7, we 
plot the output fidelity from each scheme by including the error At of beam-splitter transmis- 
sivity. Compared with Fig. 6, it turns out that a high output fidelity is still achievable and that 
the second scheme is particularly insensitive to the beam-splitter error. 
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Fig. 7. Fidelity between the ideal state Cb|0} a |0)j + C\ \ \) a \ \)b and the output state with 
the error At, = ±0.01 of the beam-splitter transmissivity = 1,2). Other parameters are 
the same as those in Fig. 6. 



6. Summary 

We have proposed two experimental schemes to generate a finite-dimensional photon number 
entangled state (PNES) with arbitrary coefficients, i.e., Y^,=QC n \n) a \n) y. One scheme is based 
on the second-order coherent superposition operation with two-mode squeezing operations, and 
the other on two first-order coherent superposition operations. We have shown that the coeffi- 
cients of the PNES can be adjusted by the parameters of beam splitters and NDPAs in each 
scheme. In particular, our schemes do not require a high-level of squeezing for the nonlinear 
materials (NDPAs) and we further demonstrated that our schemes can generate the PNESs with 
high fidelity using realistic on-off photodetectors with nonideal efficiency. The class of PNES 
is useful for CV quantum informatics as we have considered its application to quantum telepor- 
tation and nonlocality test. We have shown that the PNES of finite dimension can surpass the 
performance of the TMSS with the level of squeezing currently available in the pulsed regime. 
Furthermore, the PNES includes a broad class of non-Gaussian entangled states together with 
the TMSS (a representative Gaussian entangled state), therefore, our schemes can also be used 
for fundamental tests of quantum physics, e.g. the robustness of Gaussian versus non-Gaussian 
entanglement under noisy environments [26. 12811291 . 
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